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Non-linear dynamics of cosmic strings with non-scaling loops 

Vitaly VanchurirO 
Department of Physics, Stanford University, Stanford, CA 94305 

At early stages the dynamics of cosmic string networks is expected to be influenced by an excessive 
production of smafl loops at the scales of initial conditions Imin- To understand the late time 
behavior we propose a very simple analytical model of strings with a non-scaling population of loops. 
The complicated non-linear dynamics is described by only a single parameter TV ^ 2/(1 — C{lmin)) 
where C{1) is a correlation function of the string tangent vectors. The model predicts an appearance 
of two new length scales: the coherence length ^ ~ t/N'^ and the cross-correlation length x ~ t/N. 
At the onset of evolution A'' ~ 10 and at late times A^ is expected to grow logarithmically due to 
cosmological stretching and emission of small loops. The very late time evolution might be modifled 
further when the gravitational back-reaction scale grows larger than Imm- 

PACS numbers: 98.80.Cq 11.27.-|-d 



I. INTRODUCTION 

Cosmic strings are predicted by many models of sym- 
metry breaking phase transitions [1] and give rise to very 
distinct and detectable signatures such as gravitational 
lensing BQ, CMB non-Gaussianities 043' gravitational 
waves and ultra- high energy cosmic rays [13, lll| . 
More recently, it was also realized that cosmic super- 
strings are usually formed at the end of brane inflation 
[12 - [l4| which opens a possibility of testing the models of 
string theory in the cosmological settings. 

In the early work on cosmic strings, it was expected 
[la . [l6| that the typical length of closed loops scales 
linearly with time, but some old numerical simulations 
showed no evidence of such scaling. After more than 
twenty years of numerical studies [3, I17h201 l23l - |2a | the 
issue of string loops is still a subject of an ongoing de- 
bate. Some groups argue that the loops are much smaller 
than horizon [l9l-|22|. when others claim that the loops 
are typically produced at the near horizon scales [23l - l26| . 
Clearly the final verdict has not yet been reached. 

In addition to numerical studies a number of ver> 
teresting analytical models were proposed [27h30| ,[3^ 
some of which claim remarkable agreements with numer- 
ical simulations [30] • However, the downside of many 
models is that they often contain phenomenological pa- 
rameters, such as the mean velocity of strings, which are 
not derived from first principles. It is always left to nu- 
merical simulations to fix the unknown parameters and 
functions, which might be a dangerous path given that 
the dynamical range of the simulation is still very short. 
Moreover all of these models do not include the non- linear 
effects of loop production in a self-cosictent way. 

In this paper we will make a first step to develop a self- 
consistent dynamical model by including a back-reaction 
from non-scaling loops. The main assumptions are: 



1) non-scaling loops are predominantly produced at the 
scales of the initial correlation length ^ Imin- 

2) large scale ~ t inter-commutations do not significantly 
affect the statistical properties on smaller scales ~ Imin- 

3) in the comoving coordinates Imin remains constant, 
when all other length scales grow linearly with time. 
These assurnption are motivated by many numerical 
simulations I20I, |23|426| as well as by analytical results 

[n m m, lip? 

The paper is organized as follows. In the second section 
we derive the most important properties of the model un- 
der assumption of slowly changing power on the smallest 
scales, and in the third section we show that this power 
decays only logarithmically. The main results are dis- 
cussed in the conclusion. 



II. NON-SCALING MODEL 

One common feature of many numerical simulation 
[23l - |26| is an excessive production of very small loops 
at the sizes of the initial correlation length Imin- This 
production does not go away quickly, and even at late 
times when strings become relatively smooth the small 
loops are abundantly produced 
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At first this seems very counterintuitive given that long 
strings exhibit full scaling 



d{t) ex C(i) oc t, 



(2) 



where d is the inter-string distance and C, is the corre- 
lation length. However, one can show that the smallest 
wiggles of size Imin are very likely to form loops on a 
passage through each other if they form a cusp [30, |3a| . 
Nambu-Goto evolution of strings is usually described 
by decomposition of a position three vector into right and 
left moving waves 



'Electronic address: ' vanchurin® Stanford . edul 



x(o-, t) 



a(<T- t) + h{(T + t) 



(3) 



with condition la' 



Ib'l 



1, where prime denotes a 



derivative with respect to a. If a' and — b' curves, corre- 
sponding to the two wiggles, intersect on a unit sphere, 
then they would annihilate with a significant probabil- 
ity -30, , 35] . Therefore, if we choose two opposite- moving 
wiggles of size Imm at random then the probability p for 
them to intersect on a unit sphere can be calculated from 
the mean amplitude of wiggles 



where 



P{1) = V2-2C(0 



C{1) ^ (a'(O)a'(O) 



is a two-point correlation function and thus. 
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For the following analysis it will be useful to assume 
that all Imin size wiggles can be grouped into 



N = 



1 



P 1 — O (imin) 



(7) 



equivalence classes (we call simply directions^ ) such that 
any pair (one left moving wiggle and one right moving 
wiggle) from the same equivalence class would form a 
loop on the passage through each other and any pair 
from distinct classes would pass through without inter- 
actions. Of course the division into equivalence classes is 
not precise, but it provides us with an intuitive picture 
of a process under consideration. 

To understand the evolution of the entire network of 
strings it is convenient to introduce a new scale, we call 
the cross-correlation length^ 



X(t) = kt, 



(8) 



which describes a distance along the string on which cor- 
relations between opposite moving segment had been es- 
tablished. More precisely, x is an average distance along 
string on which the two sets: directions of all left-moving 
wave and directions of all right-moving segments, are dis- 
joint. The cross-correlation length x{t) should be smaller 
than time t due to causality, but larger than the coher- 
ence length 



m = ct, 



(9) 



which is defined as an average distance along the string 
where the direction does not change. (Note that both 
scales: coherence length and cross-correlation length are 
taken to be linear functions of time due to our original 



^ Strictly speaking each class is not a single directions, but a set 

of directions described by area ~ TrP{lmin) on a unit s phe re. 
^ The existence of cross-correlations is not new (See Ref. |27|| ). 



assumptions listed in the introduction.) If the left and 
right moving waves were completely random then one 
would have to make ~ N comparisons before an arbitrary 
pair of opposite moving directions could coincide. This 
would make x(i) ^ ^(t)-\/]V, but because of the cross- 
correlations, as we will argue below, x(i) ^ £,{t)N. 

If the decay process of a long piece of string L{t) into 
small loops is effective only on distances comparable to 
cross-correlation length then 



dLit) L{t) 



(10) 



dt ktN' 

which reflects that a given wiggle has to move a distance 
^ x(0 before meeting a "truly" random opposite mov- 
ing wiggle with whom he can annihilate with probability 
1/A^. This also tells us that an arbitrary segment has to 
move a distance ~ ^xi^) in order to form a non-scaling 
loop with a significant probability ^ 0{1). The solution 
of the differential equation (flUl) is given by 



(11) 



L ex i *=" . 
Consider a very large loop of size 






(12) 



at time t ~ ti and evolve it forward in time until time ^2 
with only non-scaling loops production turned on. Un- 
til the time when its length is larger than the cross- 
correlation length, the evolution is described on average 
by pU]) with solution ([TT|). This implies that at time 
t2 the remaining length is equal to the cross-correlation 
length 



L{t2) - kt2, 



(13) 



and afterwards t ^ t2 the decay process stops L{t) ^ 

L{t2). 

The total number of coherence segments S,{ti) — cti 
on the original piece of string is ^ — ^ 



Ncti 



-, or on average 
in each of N directions. Due to statistical fluctu- 



ations this number for left - and right-moving waves is 
not the same and under assumption of Poission statis- 
tics ^ ^vivrf '^^ ^^^ segments would survive after the 
cross-correlations are established on the entire piece of 
string L{t2). In other words 
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(15) 



Since the above result is expected to hold for an arbitrary 
ti and t2, we conclude that 



x{t) = kt 



t 

N 



(16) 



and 



or 



m = ct 

Moreover Eq. (fTOj) becomes 



N2- 



dL{t) _ L{t) 
dt ^ T' 



(17) 



(18) 



with solution L{t) oc 1/i. 

Before we move on let us quickly describe an equivalent 
model which we call the " dating model" . Consider a very 
long lineup of L men and L women moving against each 
other looking to marry someone of his (her) own type 
under assumption that there are only N distinct types 
of men and women. Once this happens a happy couple 
is removed from the lineup (they get married, have kids 
and never go back to dating). It is a simple numerical 
exercise to check that the evolution of the lineup is de- 
scribe by equations ([TO)). (IT7|) and ([T5)) . If initially men 
and women are distributed randomly, then at the very 
first step the probability to find the right person is 1/iV, 
but in the long run the probability goes down because of 
cross-correlations. So if someone did not get married by 
some late time, then he will be surrounded by "friends" 
that share the same type and gender, but when he does 
get married all of his unmarried friends are likely to get 
married as well. In the language of strings we are likely 
to see very large bursts of small loops moving with large 
velocities in the same direction as was previously argued 
in [3i, 33] . Of course the dating model does not describe 
all of the important dynamics on all scales, but neverthe- 
less it provides us with a very simple and intuitive picture 
of non-linear dynamics. All of the additional effects such 
as cosmological stretching and back-reaction of loops will 
be analyzed in the following section. 



III. EVOLUTION OF POWER 

A key assumption of the previous section was that the 
number of directions N does not change in time. Since 
this number is closely related to the average amplitude of 
wiggles P{lmin) ~ 1/VjV we must check that P{lmin.,t) 
does not vary too fast. There are two main mechanism 
which could lead to the evolution of power on the smallest 
scales: emission of loops and cosmological stretching. 

The first mechanism is caused by an enhanced produc- 
tion of small loops in the vicinity of cusps. Any left- 
moving wiggle meets N{t)'^2dt/t distinct right-moving 
wiggles per time dt, but because of the cross-correlations 
only 2N{t)dt/t ~ 2P{t)~'^dt/t of them are pointing in 
a random directions. Moreover the probability to have 
a significant overlap, in which case a significant part of 
the wiggle will be removed, is proportional to P{t)^- As 
a result, the average amplitude on scales of Imin will be 
decreasing according to 



dP(t)2 P(i)4 

ex 



dt 



(19) 



dP{t) P{tf 

ex 



dt 



t 



(20) 



In other words the wiggles with larger amplitudes are 
more likely to form cusps, and therefore loops, which 
results in an overall decrease of power on the scales com- 
parable to the loop sizes. 

The second mechanism that might be responsible for 
the reduction of power is cosmological stretching (sj 
which is described by the following equation 



P{t2) = P{tl) 



a{t2] 



2{i,^>-l 



(21) 



where a{t) ex t" is the scale factor and (w^) is 
the mean square velocity of strings. If we neglect 
the cross-correlations between opposite moving waves 
(a'(cr)b'(cr)) = 0, then (w^) would be 1/2, but as we 
now know production of small loops introduces signif- 
icant cross-correlations. Such cross-correlations can be 
roughly estimated as 



2{v^) 1 = -(a'(a)b'(c7)) ^ - 



1 



N{t) 



-P{tf (22) 



This is because the opposite moving segments can point 
in all directions but one. By combining (l?T|) and ([2^ we 
arrive to a simple differential equation which describes 
stretching of wiggles on cosmological backgrounds 



dP{t) P{tf 

ex —a- 



dt 



t 



(23) 



and has a from identical to (pOj) . 

It follows that a complete evolution equation with both 
effects (stretching and emission of loops) taken into ac- 
count is obtained by combining ([25)) and ([20]): 



dP{t) 
dt 



{A + aB) 



t 



(24) 



where A and B are some constants of order one. The 
corresponding solution is 



P{t) = {C + 2{A + aB)\og{t)y 



N{t) ^C + 2{A + aB)\og{t). 



(25) 



(26) 



This shows that the number of directions grows only log- 
arithmically which is self-consistent with our original as- 
sumptions. 

We are now ready to understand the non-scaling be- 
havior observed in many numerical simulations. At the 
initial conditions the consecutive segments of size Imin 
have large angles ~ 7r/4 which makes them an easy tar- 
get. The corresponding TV ^ 10 and the relevant length 
scales grow as 



X^O.li 



(27) 



and 



^-O.Oli 



(28) 



according to (fT6|) and (fT7|) . We would like to note that 
at early times the three length scales: correlation length 
(^, inter-string distance d and cross-correlation length x 
are approximately of the same order, but the coherence 
length is one order of magnitude smaller. Therefore, 
throughout the evolution one can neglect possible inter- 
sections of the nearby strings and production of large 
loops which would not significantly modify correlations 
and cross-correlations established on the string. In this 
limit our non-scaling model can describe fairly well the 
non-scaling evolution of wiggles and production of small 
loops. 

In general the decay of string length in a scaling net- 
work must be given by 



dL{t) 
dt 



m 



m 



t 



t 



friction ~ 



L{t) 



t 



(29) 



with Tsrnaii + ^ large + '^friction == 2. The three tcrms 
describe the transfer of energy to small loops, large loops 
and Hubble friction respectively. From ([T5)) the decay 
parameter of small loops is exactly one, but because of 
the logarithmic growth (p6)) it might be a bit smaller 
^ small ^ 1- Since the decay rate of the small scale 
power due to stretching Q is described by the mean 
velocity of strings dHJ, in the presence of very strong 
cross-correlations predicted by our non-scaling model we 
get 



Tf^cuon = -a(2(w2) - 1) ^ aP{tf < 1. 



(30) 



which is already subdominant at early times and would 
decay logarithmically as more structures on small scales 
are eliminated. Therefore our analysis suggests that in 
addition to cosmological stretching and non-scaling small 
loops we must have a population of intermediate or large 
loops with Tiarge ^ 1 which will be analyzed extensively 
in a forthcoming publication [35| . 



non-scaling loops taken into account. This was larg ely 
motivated by analytical results obtained in [30, [3l| by 
perturbative methods which cannot be used to study the 
production of loops at larger scales ^ Imin ■ The model is 
far from being complete, as it does not include all of the 
important effects on the intermediate and large scales, 
but is mature enough to make a number of interesting 
predictions. First of all the model predicts an existence 
of two new length scales: the coherence length and the 
cross-correlation length. Secondly, the model shows that 
the small scales power decays only logarithmically and as 
a result the smallest scales must play a key role in the 
dynamics of cosmic strings. 

Throughout the paper we have explicitly assumed 
that the gravitational back-reaction scale remains always 
smaller than the scale of the initial conditions Imin- How- 
ever, if the back-reaction scale can grow linearly with 
time (which may or may not be the case) then this as- 
sumption is likely to be broken sooner or later. From that 
point on, the relevant production scale of small loops Imin 
might also start to scale, but the correlation and cross- 
correlation properties analyzed in the paper might not 
be significantly modified given that Imin <C t. We expect 
that the only relevant consequence of this scenario would 
be that the primary scale of small loops emitted at cusps 
would also scale with time. On the other hand the log- 
arithmic decay of power at the scales of Imin derived in 
the previous section cannot continue indefinitely and the 
network must enter full scaling. We plan to address this 
matter in the future work. 
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